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Comparison  of  Difforont  Modes  of  Representation  of  Multivariable  Systems 

Kumpati  S.  Narendra  and  Lyle  E.  McBrido,  Jr.  . 

Gordon  McKay  Laboratory,  Harvard^University 
Cambridge,  Massachusetts 


■  'ABSTRACT  '  ■" 

A  multivariable  system  may  be  conveniently  represented  either  by  a  ma« 
trix  of  transfer  functions  or  by  a  set  of  first-order  differential  equations.  The 
mode  of  representation  depends  on  the  nature  of  the  problem  and  the  constraints 
involved,  The-^paper-deals-with  the  relative  merits  of  several  methods  as  well 
as  their  inherent  limitations^  Concepts  such  as  "controllability, "  ''observability 
'^‘structure**  and  "Interaction'*  which  are  peculiar  to  multivariable  systems  are 
examined  to  form  a  basis  for  the  comparison. 

y 


I.  Introduction  \ 

The  extension  of  single-variable  control  theory  to  linear  systems  having 
several  inputs  and  outputs  by  defining  suitable  input  and  output  vectors  and  trans 
fer  matrices  has  been  considered  by  many  authors.  Such  matrices,  defining  the 
over-all  characteristics  of  the  system,  are  logical  extensions  of  familiar  fre¬ 
quency  domain  methods  used  in  the  case  of  single-variable  systems. 

The  inadequacy  of  the  transfer  matrix  approach  to  define  Completely  even 
the  terminal  behavior  under  certain  conditions  was  pointed  out  by  Freeman[l] 
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and  Mesarovic  [2],  The  lattor  used  the  concept  of  structure  to  develop  different 
subsystem  intercohnections  to  achieve  the  same  terminal  behavior. 

study  of  optimal  control,  interest  has  shifted  to  the 
time  domain,  and  the  state  vector  method  involving  a  set  of  flrst-brder  differen¬ 
tial  equations  has  been  used  for  the  representation  of  the  rhultivariable  system. 
Terms  such  as  ''controllability''  and  "observability'*  have  been  defined  in  this 
context  to  describe  pertinent  characteristics. 

The  varipus  concepts  such  as  "structure, "  "interaction, "  "controllability, " 
"observability,  "  etc. ,  are  peculiar  to  multivariable  systems  and  the  aim  of  the 
paper  is  to  indicate  the  relation  between  these  concepts  as  well  as  the  advantages 
and  limitations  of  the  two  methods.  At  this  time,  when  the  synthesis  of  multi- 
variable  systems  is  atill  in  its  Infancy,  it  is  felt  that  a  discussion  of  the  various 
methods  of  representation  and  the  problems  to  which  each  is  best  suited  will 
facilitate  further  work  in  the  field. 

11.  Methods  of  Representation 
a.  The  Transfer  Matrix  Approach  to  Multivariable  Systems 

Over  the  last  quarter-century  a  considerable  body  of  knowledge  about 
linear  Control  systems  has  been  built  up.  Frequency  response,  transient  re¬ 
sponse  and  stability  criteria,  pole-zero,  root-locus  and  compensation  methods, 
along  with  many  other  techniques,  have  made  possible  the  analytic  design  of 
even  extremely  complex  controls  for  dynamic  systems. 

These  tools  are  based  on  the  concept  of  a  transfer  function  which  uniquely 
describes  a  real  physical  system  under  certain  conditions  of  operation.  Among 
the  conditions  which  must  be  fulfilled  if  the  transfer  function  representation  is 
to  be  accurate  are  the  following; 
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(1)  The  ayftom  ihuflt  be  approximately  linear  around  its  operating  point 
so  that  it  can  be  deicribed  by  a  linear,  iUtibhary  differeh^ial  equation  trana- 
formable  into  an  algebraic  equation  in  the  Laplace  tranaform  variable  a . 

(2)  The  ayatem  muat  be  operated  in  auch  a  manner  that  a  aignal  entere 
the  ayatem  or  a  aubayatom  at  a  preacribed  input  terminal,  and  leaves  only  at 
the  output;  information  paaaea  through  the  ayatem  in  only  one  direction,  and  a 
single  tranafer  function  deacribea  a  single ‘input,  single-output  system. 

(3)  Only  problems  in  which  all  initial  conditions  are  aero  or  of  no  im¬ 
portance  (as  in  steady-state  behavior)  can  be  treated  by  the  transfer  function 
method;  in  deriving  the  transfer  function  ail  information  about  the  effect  of 
initial  conditions  which  is  inherent  in  the  differential  equation  is  discarded. 

IThis  does  not  imply  that  the  Laplace  transform  equivalent  of  a  differential  equa¬ 
tion  contains  less  information  than  the  equation  itself;  it  is  only  in  obtatning  the 
transfer  function  that  the  terms  representing  the  effect  of  initial  conditions  are 
dropped.) 

Transfer  function  techniques  have  been  extended  to  systems  having  mul¬ 
tiple  inputs  and  outputs  by  defining  suitable  input  and  output  vectors  and  transfer 
matrices.  Such  matrices  have  been  used  in  determining  the  stability  of  muiti- 
variablc  systems,  in  designing  controllers  to  achieve  desired  dynamic  perform¬ 
ance  and  in  specifying  the  conditions  under  which  a  desired  transfer  character¬ 
istic  is  realizable. 

The  basis  of  this  extension  is  to  define  a  vector  X(s)  consisting  of  the 
Laplace  transforms  of  all  the  input  variables  in  a  given  order,  a  vector  Y(s) 
consisting  of  the  transforms  of  all  the  output  variables  similarly  ordered,  and 
a  transfer  matrix  A{s)  such  that  the  matrix  equation 
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\:  rA  systom  rdproafintad  by  a  tranitlor  matrix  A  may  be  considered  to  be  a 
combination  orslngle-varlable  systems  Interconnected  In  the  manner  of  Fig.  1. 


This  combination  of  3 ingle -variable  subsystems  is  now  a  multivariable 
system  which  obeys  Eq.  1.  If  any  combination  of  input  functions  x^(s)  is  given, 
the  resulting  outputs  y^(s)  can  be  found  from  this  equation,  since  aj^.'s  are  fixed 
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charactertatici  of  tho  aystom.  If  the  x^' a  are  the  only  yariablea  directly  affected 
by  agencies  outaide  the  ayatem,  and  the  y^'a  are  the  only  variables  bbaerved,  the 
matrix  A  contains  the  aamo  information  about  the  multivariable  syatent  as  the 
tranafor  function  doba  for  a  a  ingle -variable  ayatem. 

Howovor/  the  multivariable  ayatem  la  found  to  poSaeaa  propertiea  of  in¬ 
terest  which  do  not  exist  in  a  single-variable  ayatem.  Some  of  these  have  been 
diacuaaed  in  the  literature  under  the  names  of  interactions  •inter  relatioha  and 
Intercoupling'.': 

(1)  Input  Interaction.  Thia  quality,  which  may  be  roughly  defined  as  the 
extent  to  which  all  inputs  affect  all  outputs,  is  a  function  of  the  conventional 
transfer  matrix  A  only.  Complete  input  non-interaction  is  defined  by  BokSenbom 
and  Hood  [3]  as  tho  conversion  of  the  n-input,  n-output  multivariable  System  into 
n  single-Variabie  systems;  in  other  words,  in  a  non-interacting  system  each 
input  affects  one  and  only  one  output.  Non-interaction  is  demonstrated  in  the 
matrix  representation  by  a  diagonal  A  matrix  (or  more  generally,  one  in  which 
each  row  and  column  contains  only  one  non-zero  component). 

(2)  Output  interaction.  This  property  of  multivariable  systems  has  been 
discussed  from  various  points  of  view,  usually  with  a  desire  to  eliminate  it  as 
far  as  possible.  Output  non-interaction  is  referred  to  by  Freeman  [1]  as  *' inde¬ 
pendent  output  restoration, "  and  is  described  in  terms  ot  the  reaction  of  the  sys¬ 
tem  to  a  non-zero  initial  condition  at  output  y.,  when  all  inputs  are  identically 
zero  and  all  other  outputs  are  initially  zero.  If  under  such  conditions  the  dis¬ 
turbance  following  the  "turning  on  of  power"  in  the  system  is  confined  to  output 
y^,  the  system  is  considered  non-interacting  with  respect  to  that  output. 

Mesarovic  [2]  defines  a  similar  property  called  "interrelation"  or  "in¬ 
teraction,  "  as  the  reaction  at  other  outputs  to  an  external  disturbance  applied  to 


TR356 


.6- 

Because  of  conclltlont  fZ)  and  (3)  which  limit  the  transfer  function  to  proh'^ 
lems  involving  uhldirectiohal  signal  flow  in  the  absenco  of  initial  conditions.  Eq. 

1  does  not  contain  any  information  about  the  reaction  of  the  systern  to  either  a 
disturbance  at  the  output  or  a  non* zero  initial  condition.  To  specify  output  In¬ 
teraction,  therefore,  additional  information  not  contained  in  the  transfer  matrix 
^Ts'-neCe-ssary.  ^ 

(i)  Output  dependence.  This  term  is  intended  to  describe  the  character¬ 
istic  called  "Interrelations  based  bn  External  Changes  of  the  System"  in  reference 
[2].  It  describes  the  possibility  of  obtaining  any  Independent  set  of  desired  output 
functions  by  suitable  manipulation  of  the  input  functions.  It  n.ay  therefore  be 
Said  that  the  outputs  are  independent  if  matrix  A  is  non-singular  when  m  a  n.  If 
m  <  n  it  is  apparent  that  at  least  n  -  m  outputs  must  be  dependent  on  the  others. 

In  fact,  in  general,  if  the  rank  of  A  (the  order  of  the  highest  non-zero  determinant 
in  A)  is  equal  to  n,  the  n  outputs  can  be  independently  varied;  if  the  rank  of 
A  +  q  s  n,  then  q  outputs  can  be  expressed  as  functions  of  the  other  n  -  q.  We 
could  then  define  the  case  where  q  =  0  as  output  independence,  and  all  other  cases 
as  output  dependence  of  the  qth-order.  (Here  output  dependence  is  equivalent  to 
"infinite  intercoupling"  in  reference  (2]. ) 

The  extreme  case  of  output  independence  is  clearly  the  same  as  the  case 
of  input  non-interaction;  the  outputs  are  most  independent  (if  the  superlative  is 
permissible  here)  when  each  depends  only  on  a  single  input, 
b.  "Structure"  of  a  Multivariable  System 

The  inability  of  the  transfer  matrix  A  to  resolve  the  question  of  output 
interaction  led  the  authors  of  references  [l]  and  [2]  to  introduce  the  concept  of 
.structure. 

Although  their  points  of  view  differ  slightly,  both  Freeman  [1]  and 
Mesarovic  [2]  define  "canonical  structures"  which  are  specific  ways,  of  siib- 
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dividing  a  multivariable  •yiterh  into  single -inputi  single -output  sybsystems  de¬ 
scribed  by  transfer  fuhctlbns.  Figure  1 ,  for  oxample,  rbprbsents  the  "P-cahon- 
ical  structure"  of  reference  [2|;  various  other  combinations,  called  V-canonical, 
H  -canonical  and  non  -canonical  structures ,  are  shown  by  Mesarovlc  to  obey 
Eq.  1  equally  well,  though  the  transfer  functions  of  the  cbmpenent  single -variablo 
systems  are  no  ionger  simply  the  Components  of  the  matrix  A. 

A  proper  interpretation  of  the  structure  is  made  to  yield  some  information 
about  the  presence  or  absence  of  output  Interactioni  If  it  Is  assumed  that  all  of 
the  blocks  of  Fig.  1  are  ideal  single-variable  systems  which  are  capable  of  trans¬ 
mitting  information  in  only  one  direction,  then  a  glance  at  the  diagram  indicates 
that  no  disturbance  at  y^  can  be  transmitted  to  any  other  output.  The  P-canonlcal 
structure  is  then  said  to  be  output  non-interacting. 


Figure  2 


The  V-canonical  structure,  on  the  other  hand  (e.  g.  ,  Fig.  2)  is  one  in 
which  every  output  variable  interacts  with  every  other.  A  disturbance  (or  ini 
tial  condition)  at  will  obviously  be  transmitted  to  y2  by  way! of  the  transfer 
functions  of  ^21  and  F^.  Such  a  structure  is  described  as  ^completely  inter- 
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A  similar  use  of  structure!  or  th«  subdivision  of  multivariable  systems 
Into  slnglo'varlable'  blocks,  has  arisen  in  connection  with  the  concept  (see  Sec¬ 
tion  li  e)  of  controllability.  For  example,  kalmah  [4]  discusses  the  uncontrol¬ 
lable  system  shown  In  Fig.  3. 


Figure  3. 


If  this  diagram  Is  Interpreted  as  representing  the  state  variable  equa¬ 


tions 


5c  1  a  -X|  +  d^u 


(2) 


^2  =  "’‘2  +  **2'* 


it  is  apparent  that  d^lcj  -  djic2  =  "{^2^1  ’  *^1**2^  that  the  combination  d^Xi  i- 
d|X2  represents  an  uncontrollable  state  variable  (one  which  cannot  be  altered 
regardless  of  the  control  applied).  On  the  other  hand,  since  the  system  of  Fig. 
4  which  is  identical  in  transfer  function  has  only  one  state  variable  and  is  com¬ 
pletely  controllable,  an  apparent  association  between  structxire  and  controlla¬ 
bility  has  been  established.  Gilbert  [5]  gives  additional  examples  of  structures 
which  are  not  completely  controllable;  he  also  emphasizes  the  fact  that  the 
transfer  matrix  does  not  necessarily  indicate  the  correct  Order  of  the  system 
by  drawing  structures  of  different  orders  which  have  the  same  transfer  matrix. 
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Figuro  4 


t.Thd  State  Vac  tor  Method 

The  frequency  domain  approach  using  the  concept  of  transfer  function^ 
described  in  the  preceding  sections/  has  been  used  in  the  design  of  control  sys¬ 
tems  for  over  two  decades.  In  spite  of  its  wide  application  the  method  suffers 
from  several  inherent  limitations.  The  tenuous  relationship  existing  between 
frequency  and  transient  responses  of  a  system  results  in  the  designer  having 
only  qualitative  information  about  transient  response.  Hence,  when  the  tran¬ 
sient  response  Is  rigidly  specified  the  frequency  domain  method  is  not  usually 
suitable.  The  limitations  of  the  method  also  become  obvious  when  dealing  with 
non-linear  or  time<varying  systems.  When  time  responses  are  of  interest,  it 
is  found  that  a  direct  investigation  of  the  behavior  of  the  differential  equations 
governing  the  system  directly  yields  more  pertinent  results.  This  has  resulted 
in  the  '*state  vector  approach"  wherein  the  differential  equations  describihg  the 
physical  system  are  studied  and  the  system  behavior  is  controlled  directly  in 
the  time  domain. 

In  the  state  vector  method  the  multivariable  system  is  represented  by  a 
set  of  differential  equations  of  the  form 

*,  =  f.(Xi.  ...  x„,  u,.  Uj  ...  u,,  Pj.  .  .  p^.  tj  i  =  1,  ....„  (3) 


xj(0)  =  C 


or  in  vector  notation 
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X  »  PtX;  Vi  P;  t)  X(0)  ■  C  / 

Tho  Xj  reproHont  "n"  state  variablos  which  c6mplotbi|^  specify  the 
plant  at  any  instant, U^(t)  are  the  reference  inputs  to  the  bys tern  and  P^(t)  the 
unintentional  dlstuabances  entering  the  system.  The  coordinate  system  with 
Xp  ^0  coordinates  spans  an  n-spaco  which  is  caliod  the  state  space. 

The  vector  X  whose  components  are  x^v  the  state  variablos.  is  called  the  state 
vector  and  tho  curve  traced  by  the  state  vector  with  tho  passage  of  time  is 
known  as  a  trajectory.  It  is  worthwhile  remembering  that  the  mathematical 
model  (represented  by  Eq.  3)  is  only  an  assumption.  The  slrripliflcatioh  intro¬ 
duced  by  the  as  is  com- 

In  general. the  problcmiO^C0mr^^^^^liT^^^M^^m^y'%e  coris'l^or'ed  as  the 


lis  approachfmayioe  considered  as  the 


problem  of  determining  the  inputs  U  tb  the  systeni  such  that  the  output  X(t)  cor¬ 
responds  as  closely  as  possible  to  some  desired  behavior. 

If  we  confine  our  attention  to  the  problem  of  control  of  systems  which 
are  linear,  the  differential  equation  of  the  system  may  be  written  as 

X(t)  =  A(t)  X(t)  +B(t)  U(t)  +F(t)  (4) 

where  X  andFare  n  dimensional  vector  functions,  A  is  an  (n  x  n)  matrix 
function  and  B  an  (n  x  r)  matrix  function.  The  very  first  question  that  arises 
in  the  analysis  of  such  a  system  is  whether  or  not  the  state  of  the  system  can 
be  controlled  arbitrarily  by  the  application  of  a  suitable  control  function  U(t), 
i.e.  ,  whether  any  state  of  the  system  can  be  driven  to  any  other  desired  state, 
by  suitable  control  action.  Such  a  system  has  been  defined  to  be  controllable 
[4].  It  is  obvious  that  in  general  controllability  is  a  prerequisite  to  any  type 
’:of  optimal 'Control.^ '■  -  ■  ' 

The  solution  of  Eq.  4  represents  the  behavior  of  the  state  of  the  system, 
and  is  defined  by  a  set  of  functions  ^j(f)-  Thus  X{t)  implies 


(5) 


X(t;  x^.  g 


where 


Xo.  g  -  x^ 


and  Xj(t)  are  the  components  of  X(t)  in  the  state  space.  If  the  functions  A(t)»  D(t) 
and  U|t)  are  defined  for  ail  *00  <  t  <  a>  and  are  bounded  for  each  t,  the  solution 


X(t)  .  i(t,  g  X^  4  J  •(t,  t)  B(t)  U{t)  dT>  ^  f(t,  t)  F<t)  di 


(6) 


where  i(t,  g  is  the  principal  matrix  solution  of 


ft iiSftf  xttf 

and  is  called  the  transition  matrix  of  Eq.  3.  In  the  absence  of  external  disturb¬ 
ances  F(t),  Eq.  4  reduces  to 

t’ 


Xft)  =  #(t,  g  X^  +  J  »(t.  t)  B(t)  U(t>  dr 

■■  ft>  ■ 


<7) 


and  when  u(t)  is  a  scalar  to 


X(t)  =  »(t,  g  »(t,  t)  B{t)  u(t)  dr 


18) 


where  B(t)  is  an  n  dimensional  vector  function.  A  state  (X  )  at  time  t  may  be 


considered  to  be  controllable  if  there  exists  a  control  U{t)  which  transfers  the 
system  to  the  state  X  =  0  at  some  time  tj 


X{t^)  =  0  =  g  4^  4(tj.  t)  B{t)  U(t)  dr 


(9) 
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io  that  -  t)  B{t)  U|t)  dr  holdi  for  ao^m  U(t).  The  mathematical 

theory  of  the  controllability  of  linear  dynamical  ayatoma  la  dlacuaaed  In  conald- 
erable  detail  In  (6)  and  the  condltlona  for  eontrollablllty  of  time  -varying  and 
^^rcohstant' aystema 'are 'de'rived. 

Conalderlng  the  caae  when  A  and  B  In  Eq.  4  are  conatant,  the  ayatem  la 
completely  controllable  If.  and  only  If,  the  conatant  n  x  rn  matrix  ^,  AB,  A^B, 

. . .  a”  “  ^  haa  a  rank  **n.  •*  Thla  ehaurea  that  the  control  function  U(t)  operated 
on  by  B  apana  the  entire  atate  apace,  ao  that  each  atate  variable  may  be  Inde- 
'^^^dentl  y^a  ffe  ct  etf^y Jhe  ^,on¥fo  ihaiifolh>^ 
neceaaary  and  aufflclent  condltlona  of  controllability  aa  well  aa  other  mathemat¬ 
ical  concepta  the  reader  la  referred  to  reference  (6). ) 

Conalderlng  the  time -Invariant  caae  the  equation  deacrlblng  the  ayatem 
la  of  the  form 

X  -  AX  +  BU  (10) 

when  U  la  a  vector,  and 

^  ^  ^(11) 
when  U  is  a  scalar  quantity.  Since  by  definition  of  controllability  it  is  clear 
that  it  is  invariant  to  any  non-singular  linear  transformation,  if  we  define 
=’PY 

.■'Eq. '10  becomes.'- 

'  -Y  ■='P"V’AP  Y.+ P’^-'B  U.  -y:  .r' V 

By  the  proper  choice  of  P  it  is  seen  that  the  equation  may  be  reduced  to  normal 
form.  If  in  Eq  10  it  is  assumed  that  X  is  already  in  its  principal  coordinates. 
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A  li  In  its  Jordan  canonical  form.  Assuming  that  the  matrix  A  has  "q'*  distinct 
eigenvalues  Xj,  ...  and  an  eigenvalue  X^  ^  ^  of  niiultipilclty  (n  -  q)i  Xf  B 
and  A  in  Eq.  11  may  bo  expressed  in  the  following  form 


X  ■  <  X. 


‘q  +  1.  1 


’q  +  1,  I 


0  I 


jV'V 

i“\  \ 


q  +  1,  n  -  qj 


q  +  1,  n  - 


The  transition  matrix  is  of  the  form 


4  (t)  = 


‘xX_t  I 

A  !  .  L  .  .  .  .  .  . 

ffC  t  tf 

e  I  0  1 


(nTinT! 


For  complete  controllability,  4  (t)B  should  have  "n”  components  which  are  line 
arly  independent  and  non-zero  over  each  positive  interval  of  time. 

The  block  diagram  representation  of  Eq.  11  is  shown  in  Fig.  5.  This 
in  turn  may  be  expressed  in  block  diagram  form  as  shown  in  Fig.  6. 


Figure  5 


To  make  $(t)Bhave  "n”  linearly  independent  components  it  is  clear  that  b^, 
b,  .  .  .  b  and  b  .  ,  ^  ^  +  1  n  -  d  ^  ® 

2  q  q  +  I,  n  -  q  HT  “  4  ; 

ensures  that  the  state  variables  ^  ^  T  •  ^q  +  1,  n  -  q  controlled 

independently.  In  other  words,  at  least  one  of  the  control  inputs  should  affect 
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X|,  . . .  and  ’‘q  ^  j  .  axpreggod  by  the  fact  that  tha  matrix 

£b,  AB,  ...  a”  *  ahoutd  havo  a  rank  **n.  ** 

ControUabllity  as  dofinad  abova  impllos  that  by  tha  aultabla  aippllcatioh 
of  a  control  function  U  the  state  variables  may  be  Affected  Independehtly  In 
time  so  that  for  some  value  of  **1**  the  variables  x^  may  be  made  to  assume  a 
speelfted  value.  From  the  above  discussion  it  is  clear  that  a  system  of  the 
form  shown  in  Fig.  7a  is  controllable  while  the  one  indicated  in  Fig.  7b  is  not. 


Figure  7 

When  the  state  vector  representation  is  used  for  the  description  of  a 
multivariable  system,  the  problem  generally  becomes  the  determination  of 
U(t)  (subject  to  constraints  on  U  and  X)  to  minimize  or  maximize  a  perform¬ 
ance  criterion.  The  state  vector  X  is  driven  from  an  initial  state  X.  to  a 
final  fixed  state  X^  or  more  generally  a  varying  state  Z(t). 

A  second  concept  [4]  which  is  of  importa.nce  while  using  the  state  vector 
method  is  termed  "observability"  and  applies  to  the  situation  when  the  state 
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vectori  are  not  directly  accesalbla  but  have  to  be  computed  from  the  available 
data,  in  particular, in  the  linear  caee, If  the  output  Y  la  such  that  Its  components 


where  C  is  a  p  x  h  matrix  function  and  Y  Is  a  p  x  I  vector.  The  system  Is 
'’observable'’  if  all  the  state  variables  Xj  of  the  staid  vector  oan  be  de^^^ 
by  observing  Y  over  a  finite  Interval  of  time. 

Assumlnjt  that  the  equations  of  the  system  are  expressed  In  terms  of  the 


principal  coordinates  and 


AX  +  BU 


X  the  block  diagram  representa 


and  if  A  has  n  distinct  eigenvalues  X^ 
tlon  of  the  system  Is  as  shown  In  Fig.  8. 
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It  la  seen  that  each  component  V  ^  pf  the  ybctor  V  controls  the  corresponding 
state  X|.  irtho  ith  row  of  B  contains  only  aero  blcmehts,  Is  aero  and  x^. 

Is  uhcontrbllable.  Similarly,  If  the  elements  pf  the  jth  column  of  C  vanish,  Xj 
Is  not  dbaorvablo-  This  definition  is  suggested  In  rofo renCe  (5). 

For  the  more  complex  case  when  the  system  has  multiple  characteris¬ 
tic  roots  SuCh  a  simple  definition  Is  not  valid.  Due  tb  the  Internal  feedback  In 
the  systern,  Fig.  5  variations  in  Xj  may  be  observed  by  bbservlng  the  bther 
state  variables.  Similarly,  the  system  may  be  completely  controllable  even 
though  a  row  of  the  B  matrix  vanishes. 

111.  Relative  Merita  of  Various  Methods  of  Representation 
.e, '"'The  Transfer  Matrix  V.  .'v/ 

The  matrix  A,  which  relates  Input  and  output  variables  In  the  frequency 
domain,  Is  the  simplest  and  most  direct  moans  of  applying  techniques  developed 
for  single -variable  servomechanisms  to  the  field  of  multivariable  systems. 

(1)  Input  non-interaction.  In  addition  to  the  problems  of  controller 
design  involving  stability,  steady-state  errors  and  similar  properties  of  the 
component  subsystems  (a  humber  of  which  have  been  studied  in  the  literature 
by  means  of  the  transfer  matrix),  the  synthesis  of  a  controller  to  produce  in¬ 
put  non-interaction  from  an  interacting  plant  may  be  highly  desirable  in  order 
to  reduce  the  number  of  degrees  of  freedom  of  the  designer  and  simplify  the 
dynamic  analysis.  The  requirements  for  non-interaction  can  be  directly  and 
generally  obtained  from  the  matrix  A. 

Let  us  first  consider  a  given  interacting  plant  of  known  transfer  ma¬ 
trix  Ap,  which  is  to  be  connected  in  cascade  (Fig.  9)  with  a  forward  acting 

controller  whose  transfer  matrix  is  A  . 

■■  c  - 
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Figure  9 


Then  A  A^X  >  Y  ■  AX.  and  A»  the  roaulling  over -all  transfer  function, 
Is  to  be  diagonal.  The  control  elements  are  determined  from 


hence 


^ij  ^  *jj  ®pij 


tsii  ,  M 

Vj 


The  control  elements  a^,,  can  be  determined  If  A  is  non-singular. 

clj  P 

From  Eqs.  17  and  18  it  is  apparent  that  the  ratios  of  the  elements  in  a  given 


column  of  the  matrix  A^  are  fixed  by  the  condition  that  the  total  system  be  non* 
interacting  (that  A  be  diagonal).  However,  the  freedom  to  determine  an  arbi¬ 
trary  transfer  matrix  A  remains,  as  each  column  may  be  multiplied  by  an 
arbitrary  transfer  function  a... 

If  control  elements  are  all  connected  in  the  feedback  path,  the  block 
diagram  of  Fig.  10  results,  and  .  l  ^ 

Ap(X  +  A^Y)  =  Y  or  =  AX 


A  =  (I  -  A  A  )A  A  -  A  =  A  A  A 

p  P  'C'-  .  .  ^p  ^  'P  .C^  •. 


A  =  A~^  -  A"^ 

'  -:C 


A'^  =  (-L  5..] 

-ij 


rrs^.  ' 
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Then  we  find  that 


In  this  case,  all  off-diagonal  elements  of  the  controller  are  fixed  by  the  non¬ 
interaction  requirement,  but  the  diagonal  terms  remain  free  to  establish  any 
desired  diagonal  transfer  matrix  A.  Here  both  A  and  A^  must  be  non-singular 
requirements  of  non-singularity  are  not  surprising.  If  Ap  were 
singular  it  would  indicate  that  one  or  more  of  the  outputs  is  not  independent 
(is  a  constant  function  of  the  other  Outputs,  regardless  of  the  combination  of 
inputs  applied).  It  is  clearly  impossible,  under  these  conditions,  to  make  the 
outputs  vary  independently  by  changing  the  plant  inputs.  If  A  Were  singular, 
it  would  indicate  that  one  output  is  always  zero  (since  A  is  diagonal).  This  is 
possible  with  a  feed-forward  controller,  since  by  a  proper  selection  of  the  n 


inputs  any  number  of  outputs  can  be  made  zero  (if  A  is  non-singular).  But  the 


feedback  controller  would  then  only  have  n  -  l  hon-zero  inputs,  and  consequently 
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ita  output  could  contain  only  n  •  1  independent  functions  of  the  plant  output,  not 
enough  to  determine  uniquely  the  system  outputs. 

Thus,  in  ah  h  x  h  system  there  always  remain  n  conditions,  a:fter  non. 
interaction  has  been  attained  which  can  be  chosen  to  produce  the  desireci  dlago. 
nai  transfer  matrix.  The  use  of  a  feedback  controller  has  the  property  (Which 
may  be  advantageous)  that  the  controller  cbmponents  determined  by  the  cohdi. 
tlon  of  nbn.interactlon  are  Independent  of  the  dynamic  characteristics  chosen 
for  the  diagonal  elements  of  A. 

(b)  The  limitations  of  5structure.  *  The  use  of  structure  as  a  necessary 
property  in  the  definition  of  a  multivariable  system  was  first  suggested  in  order 
to  provide  information  about  output  interaction:  if  a  multivariable  Siystem  is 
really  made  up  of  ideal,  isolated,  unidirectional  single^arlable  systems  (as 
the  subdivided  block  diagram  implies),  then  it  is  clear  that  the  system  of  Fig. 

1  is  output  non-interacting  while  that  of  Fig.  2  is  Completely  interacting. 

A  difficulty  may  arise,  however,  if  we  require  quantitative  information 
about  the  Interaction.  If  all  Inputs  are  identically  zero  (X  =  '  0)  and  a  known 
disturbance  y|(8)  is  applied,  what  is  the  resulting  y^i^)?  in  the  system  of  Fig. 
1,  the  answer  is  obviously  y2  S0.  In  the  structure  of  Fig.  11 
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72  *  Vi  V -canonical  atructure  of  Fig.  2»  while  It  aeami 

clear  that  y2(*)  /  0,  the  only  equatlona  at  our  cllipoial  (when  X  ■  0)  are 

0  •>  V,jFjy,  *y^  . 

or  In  matrix  form  (1  -  FV)Y  ■  0,  which  in  general  has  only  the  trivial  solution 

The  problem  Is  directly  traceable  to  condUton  (2)  which  was  stated  as 
a  requirement  for  transfer  function  representation  In  Section  II:  signal  flow 
must  be  unidirectional.  We  have  been  able  to  violate  this  condition  In  the  case 


of  Figs.  1  and  1 1  only  by  tacitly  defining  the  reverse  transfer  function  Of  a 
block  to  be  zero.  As  long  as  a  given  block  encounters  Input  signals  at  only  one 
of  Its  terminals,  this  assumption  provides  a  solution;  In  Fig.  2,  however, 
block  F|  encounters  an  Input  signal  simultaneously  at  both  terminals:  y^  at 
the  output  and  y[V21^2^12  input.  Clearly,  no  solution  Is  possible  In 

the  general  case.  In  order  to  obtain  a  meaningful  quantitative  description  of 
output  interaction,  a  system  description  valid  for  bidirectional  signal  flow  Is 
needed.  Such  a  representation  is  discussed  in  the  following  section. 

A  similar  difficulty  arises  in  the  representation  of  uncontrollable  (or 
unobservable)  systems  by  means  of  a  structure,  or  block  diagram  of  first- 
order  transfer  functions.  As  mentioned  in  Section  II,  an  intuitive  interpreta¬ 
tion  of  Fig.  3  results  in  the  equation 


where  z  =  *12^1 


■*1*2- 


Consider,  however,  the  network  of  Fig.  12  where  u  is  the  input  voltage 


and  and  x^  are  output  voltages  measured  at  the  points  indicated. 
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Figure  12 


If  the  component  valuee  are  property  choeen/  the  applicable  equatione  are 


where  a  » 


X2«  If  these  equations  are  transformed  into  the  frequency  domain* 


and  all  initial  conditions  arc  set  equal  to  zero,  the  resulting  transfer  function 


equations  are 


corresponding  to  the  block  diagram  of  Fig.  13*  which  is  similar  to  that  of  Fig.  3 


Figure  13 


1 

s  1 
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In  this  cat«,  however,  th« re  are  two  ihittal  condlttons,  3i(0)  and  fc(0),  which 
cannot  be  brought  to  zero  by  any  choice  of  u. 

From  this  example  it  can  be  seen  that,  becauso  the  transfer  function 
is  valid  only  when  all  initial  conditions  are  zero,  neither  a  transfer  function 
nor  a  bibck  diagram  made  up  of  transfer  functidna  can  uniquely  specify  the 
number  of  ihltlal  eonditions  necessary  to  define  a  general  ablutlbn  to  tite  sys¬ 
tem  equations,  and  therefore  neither  the  transfer  matrix  nor  the  structure 
can  specify  the  order  of  the  system,  although  the  minimum  order  can  be  com¬ 
puted  from  the  transfer  matrix.  The  statement  of  reference  [5]  that  uncontrol¬ 
lable  or  unobservable  elements  do  not  appear  in  the  transfer  matrix,  applies 
equally  to  the  block  diagram  or  structure  of  the  system. 

In  short,  the  Structure  or  subdivision  into  single-variable  blocks  of  a 


multivariable  system  gives  no  more  information  about  the  response  of  the  sys¬ 
tem  than  docs  the  transfer  matrix.  This  result  is  well  known. in  the  case  of 
single<wariable  systems.  • 

The  state  vector,  Or  differential  equation  representation,  which  de¬ 
scribes  completely  the  effect  of  initial  conditions,  can  be  used  to  determine 
whether  or  not  output  interaction  exists.  Since  1^  the  response  of  x.  to 

a  unit  initial  condition  on  x^,  Complete  output  non-interaction  in  accordance 
with  Freeman's  definition  corresponds  to  a  diagonal  matrix;  if  (j>.j(t)  ~  0 


an  initial  condition  in  has  no  effect  on  x^(t). 

While  output  Interaction  of  this  type  can  be  expressed  in  terms  of  the 
4  matrix,  if  interaction  is  more  generally  interpreted  as  the  effect  of  an 
externally  imposed  change  in  one  output  variable  on  another  output  the  state 
vector  representation  does  not  provide  adequate  information. 
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c.  A  Multivariable  Syatom  aa  a  Gonoralizod  Electric  Network 


One  multivariable  ayetem  which  hat  been  atucliod  In  great  detail  ia  the 
electric  network  made  up  of  conatant  impodancea  and  atatiohary  voltage  aourcba. 
In  additlont  by  means  of  either  the  Impedance  analog  or  the  mobility  analogi 
mechanical  systems  can  be  expressed  In  terms  of  equivalent  hetworka/  F^^^ 
more)  e  great  many,  if  not  allt  other  types  of  linear,  atatlonary  aystems  are 
governed  by  the  same  differential  equations  applicable  to  networks ;  thus  a  wide 
variety  of  systems,  including  all  linear  electromechanical  systems,  can  be 
studied  by  means  of  network  analysis. 

(1)  Inputs  and  outputs.  To  make  a  connection  between  the  representa¬ 
tion  of  the  preceding  section  and  the  network  representation,  let  us  consider  a 
system  which  has  n  inputs  X|  .  . .  x^  and  moutputs  y^  •  •  •  Vu,* 

In  the  above  discussion  it  was  assumed  that  the  inputs  x^  are  independent 
variables  which  may  be  manipulated  for  the  purpose  of  altering  the  outputs  y^; 
the  outputs  are  defined  as  variables  which  cannot  be  directly  manipulated  ex¬ 
ternally,  but  in  whose  values  we  are  interested.  It  is  assumed  that  although 
other  variables  exist  inside  the  system  they  cannot  be  manipulated  or  are  not 
ofinterest;  therefore,  they  are  omitted. 

If  the  system  is  an  electromechanical  system  with  n  terminals  (or  more 
properly  h  ports)  it  is  apparent  that  two  variables  are  present  at  each  terminal. 
For  example,  a  single  shaft  is  associated  simultaneously  with  a  position  and  a 
torque;  an  electrical  terminal  similarly  has  a  voltage  and  a  current.  Other 
combinations  of  variables  may  be  chosen  (velocity  and  torque,  voltage  and  pow¬ 
er,  etc.  )  but  the  condition  remains  that  the  state  of  one  terminal  is  specified  by 
two  variables.  Thus  an  n-port  network  has  2n  variables  existing  at  its  terminals 
If,  for  example,  the  n  voltages  are  fixed  by  external  connections,  the  n 
currents  are  completely  determined  for  a  given  network.  In  general,  no  current 
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can  be  altered  without  simultanoousiy  changing  one  or  more  voltages.  We  can 
then  say  that  the  n  voltages  are  Independent  variablesV  Tneaning  that  they  are 
dotormlhod  by  force*  outsidei  the  hetwork;  the  n  currents  must  then  be  dcpend- 
ent»  determined  by  the  network  charactorlstlcs  and  by  the  n  independeht  volt¬ 
ages,  Of  course  it  is  not  necessary  that  the  independent  variabies  all  be  volt¬ 
ages  dr  all  currehtsi  any  combination  of  voltages  and  currehtsmay  be  consid¬ 
ered  Independent^  with  one  restrictidhs  the  voltage  and  Current  at  any  one  ter¬ 
minal  may  not  simultaneously  be  both  independent  or  both  dependent.  This  is 
a  consequence  of  the  general  statement  that  voltage  and  current  at  any  point  are 
determined  by  the  simultaneous  solution  of  two  equations,  one  representing  the 
current-voltage  relationship  on  one  side,  the  Other  representing  a  similar  re¬ 
lation  on  the  other  side  of  the  point.  Neither  conditions  Inside  the  network  nor  ; 
those  outside  can  alone  determine  both  current  and  voltage. 

If  we nowconsider  this  network  as  a  control  system,  we  must  classify 
the  inputs  and  outputs  as  dependent  or  independent.  Clearly  the  inputs  are  in¬ 
dependent,  since  they  are  manipulated  from  the  outside,  and  the  outputs  are 
dependent,  since  they  are  to  be  Controlled  by  the  inputs.  Since,  however,  the 
inputs  are  not  generally  found  at  the  same  terminals  as  the  outputs.^he  selec¬ 
tion  of  n  inputs  and  m  outputs  will  mean  that  the  system  has  n  i  m  ports,  and 
hence  2(n  -f  m)  terminal  variables.  We  then  are  left  with  m  independent  vari¬ 
ables  which  are  not  inputs,  and  n  dependent  variables  which  are  not  outputs. 
These  correspond  to  the  conjugate  variables  at  the  output  and  input  terminals. 

Let  us  call  the  ni  independent  variables  secondary  inputs,  and  the  n  dependent 
variables  secondary  outputs. 

(2)  The  Kron  impedance  mixed  method.  A  very  general  method  of  solu¬ 
tion  of  electric  network  problems  is  the  Kron  mixed  method  described  in  Chapter 
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IV  of  reference  [7].  It  is  applicable  in  its  original  form  to  networks  containing 
any  combinatlbn  brconstant  inductances,  constant  capacltanCeSt  constant  resist' 
ances,  and  constant  frequency  voltage  and  current  generators.  As  external 
perturbations  or  Inputs/it  is  assumed  that  between  any  pair  of  ternalnale  there 
m.’iy  be  connected  eitho r  a  constant  voltage  source  or  a  constant  current  eource . 
It  Is  assumed  that  if  the  number  of  external  voltage  and  curront  sources  Is  less 
than  the  humber  of  node  pairs  of  the  network,  the  remaining  node  pairs  are  con¬ 
nected  by  voitmeters.  A  "voltmeter  mesh"  may  be  considered  equivalent  to  a 
constant  current  generator  of  aero  amplitude. 

If  our  network  is  to  represent  a  system  which  has  a  transfer  matrix  of 
the  type  discussed  In  the  previous  section,  all  outputs  must  be  xero  wh^i  all 
Inputs  (primary  and  Secondary)  are  zero.  For  the  network,  this  condition  is 
true  only  if  no  internal  sources  are  present.  Our  basic  network  thus  Contains 
only  impedances.  This  condition  docs  not  exclude  mutual  impedance  between 
branches,  whether  symmetric,  asymmetric  or  unilateral. 

The  network  impedances  are  defined  in  terms  of  the  relation  e  s  a(u)l» 
where  e  and  1  are  the  complex  descriptions  of  a  sinusoidal  voltage  and  current 
of  frequency  w.  This,  however,  is  commonly  recognized  as  being  a  special 
case  of  c(s)  =  z{s)i(s)  where  the  transform  variable  s  is  replaced  by  ju.  We 
can,  therefore,  interpret  the  impedances  as  transfer  functions,  and  after  solv¬ 
ing  the  network  equations  find  the  output  currents,  for  example,  for  any  arbi¬ 
trary  input  form. 

The  notation  in  the  following  discussion  is  in  part  that  of  reference  [7], 

If  there  are  no  internal  sources,  the  network  obeys  the  matrix  equation 
E'  =  Z'j’  (from  reference  [7],  Eq.  44-10) 


(24) 


TR356 


-27- 


"Ol  ■ 

J'- 

■*Jl  ■ 

»  ■  ■  ■  ■ 

k. 

• 

Ey 

•  • 

V 1 1 
• 

J 

2'-  CjZC 


(25) 


2  is  the  network  impedance  matrix,  C  is  the  transformation  matrix  re¬ 
lating  the  branch  currents  to  the  mesh  currents,  is  the  transpose  of  C,  is 
the  known  voltage  applied  to  the  ith  node  pair,  *Vj  is  the  unknown  voltage  appear 
ing  across  the  jth  node  pair,  the  unknown  current  in  the  kth  Maxwell  mesh 

*i^  is  the  unknown  current  in  the  external  voltage  generator,  E^  and  Ij  is  the 
known  current  applied  in  the  external  circuit  between  the.  nodes  of  the  jth  pair. 
Also  the  subscript  Mis  the  number  of  Maxwell  meshes  in  the  original  netwdrk, 

V  is  the  number  of  external  voltage  sources  applied  and  P  is  the  number  oX  node 
pairs  of  the  network.  (V  <  P  is  required  by  Kirchoff's  laws.) 

In  the  above  equations,  the  starred  scalar  quantities  are  unknown,  the 
matrices  Z  and  C,  and  hence  Z',  are  characteristics  of  the  network,  the  applied 
voltasjes  and  currents  E.  and  I.  are  known  and  the  first  M  elements  of  E'  are 


zero  because  of  the  absence  of  internal  voltage  sources.  ^ 

We  can  now  assume  that  our  inputs  are  among  the  applied  external  vari¬ 
ables  E^  and  I.,  and  that  our  outputs  are  among  the  starred  unknowns.  Since  Eq. 
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24  ropresonts  P  >  M  equations  in  P  +  M  unknowns ^  our  outputs  can  be  solved 
for  In  terms  of  the  Inputs  and  other  externally  determined  variables,  but  we  do 
not  have  an  explicit  transfer  matrix.  Also,  Eq.  24  Invoivos  a  number  of  mesh 
currents  which  are  completely  inte  rnal  and  therefore  of  no  Interest  to  us  in 
describing  the  transfer  characteristic  of  the  network;  there  may  also  be  a  num< 
ber  of  node  pairs  where  no  oxternal  connection  has  been  made  and  the  vbltag;es 
across  these  are  of  no  Interest. 

In  writing  the  vectors  E'  and  J'  above,  the  voltages  were  arbitrarily 
ordered  so  that  known  variables  appear  In  a  block,  and  the  unknowns  In  an 
adjacent  block;  however,  Eq.  24  holds  regardless  of  which  voltages  are  known. 
As  long  as  P  voltages  and  currents  are  given  (any  number  may  be  Identically 
sero)  the  P  i*  M  unknown  voltages  and  currents  may  be  found. 

By  partitioning  the  Z*  matrix  and  the  voltage  and  current  vectors,  and 
performing  the  necessary  matrix  operations  (Involving  inversions  which  are 
justified  in  reference  [7]),  an  equation  of  the  following  form  can  be  obtained: 


where  Q  la  a  P  X  P  matrix  derived  from  Z'  by  partition,  appropriate  matrix 
operations  bn  the  submatrices,  and  reassembly.  v 

Equation  26  represents  the  complete  terminal  solution  for  this  partic¬ 
ular  choice  of  dependent  and  independent  variables.  If  this  choice  is  altered, 
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the  manner  in  which  Z*  if  partitioned  change!  and  therefore  the  reiultant  Q  if 

'^.changed.".".'  " 

If  no  connection  if  to  be  made  to  a  certain  node  pair  (terminal)  the  cor - 
refpOnding  current  If  iet  equal  to  aero,  and  the  correapohdlhg  voltage  if  omitted 
af  being  of  ho  intereft.  Thif  ellminatef  one  row  and  ohe  column  from  the  ma- 
trik  Q.  Thif  can  be  continued  until  only  the  terminalf  In  which  we  are  Interefted 
remain.  We  now  have  a  matrix  equation.of  the  form: 

D*  cy  I  ^  .  .. 

where  D  if  the  vector  of  all  dependent  variahlef  (primary  and  fecondary  outputf), 
and  Ilf  the  vector  of  all  independent  variablea  (primary  and  fecondary  Inputf) 
and  Q'  if  the  matrix  Q  with  all  rowf  and  columnf  not  correaponding  to  input  or 
output  terminalf  omitted. 

We  have  now  obtained  a  transfer  function  involving  ^  terminal  variables, 
for  a  linear  electric  network.  If  the  equations  governing  a  linear  multivariable 
system  are  put  into  this  form,  will  thif  give  us  additional  information  about  its 
behavior?  Note  that  this  form  can  be  obtained  only  after  it  has  been  decided 
which  variables  are  independent  and  which  are  dependent.  The  transfer  matrix, 
like  the  transfer  function,  remains  unilateral  in  the  sense  that  the  role  of  in¬ 
dependent  and  dependent  variables  cannot  be  Interchanged  without  altering  the 
transfer  characteristic.  It  is,  however,  bilateral  in  the  sense  that  a  given 
terminal  can  simultaneously  receive  a  signal  and  transmit  one. 
d.  A  Transfer  Matrix  Representation  of  a  Multivariable  System 

Let  US  define  the  vectors.  Here  V.' 3  are  the  dependent  variables,  u.'s 
are  independent  variables,  y^'s  are  true  outputs,  Xj^'s  are  true  inputs,  is  the 
secondary  output  conjugate  to  (existing  at  the  same  terminal  as)  x^,  and  w.  is 
the  secondary  input  conjugate  to  y.. 


TR356 


-30- 


Now  an  equation  of  the  eame  form  aa  27,  but  deacrlblng  an  n  -f  m  terminal 
multivariable  ayatem  can  be  written  at  ‘  , 


^{m  +  1)1 


^(m  +  n)  1 


^{m  +  l)n 


^(m  +  n)n 


(m  +  l)(n+  1) 


(m  +n)(n+  1) 


'  ^(m  +  l)(m  +n) 


’  f(m  +n)(m  +  n) 
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Then  Eq.  29  is  equivalent  to  the  two  equations  , 

Y  =  AX  +  BW  (30) 

■jinii'''';: 

, ■  ■  Z  M  CX..+  DW 

If  wo  assume  that  the  Independent  variables  w^  conjugate  to  the  outputs 
are  all  equal  to  aerot  then  Eq.  30  Is  rbduced  to  Eq[.  1;  But  this  is  JuVt  w^ 
is  ordinarily  done  in  single -variable  theory;  if  the  output  is  a  position,  all 
known  loads  attached  to  the  output  shaft  are  considered  to  be  a  part  of  the  sys¬ 
tem,  and  the  externally  applied  torque  at  the  outptrt  shaft  is  therefore  assumed 
to  be  aero.  The  two  independent  variables  are  thereby  reduced  to  one,  the  in¬ 
put.  Such  a  system  is  therefore  a  special  case  of  a  multivariable  system  which 
has  its  secondary  input  equal  to  aero. 

If  the  network  described  by  the  matrix  P  is  the  entire  system  under  con¬ 
sideration,  the  secondary  outputs  a^  are  of  no  Interest,  and  Eq.  31  can  be 
ignored.  If,  however,  we  wish  to  interconnect  this  network  with  another  de¬ 
scribed  by  the  matrix  P',a|  acquires  a  hew  importance.  Suppose  the  input 
signal  X|  is  to  be  obtained  from  the  output  y'^  of  network  PV.  If  and  y^  are 
both  voltages,  this  result  can  be  obtained  by  connecting  terminal  1  to  terminal 
1*.  This  connection  requires  at  the  same  time  that  the  two  currents  be  equal, 
or  that  a.  =  w'. .  The  corresponding  block  diagram  is: 


Figure  14 
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Now,  even  though  only  y|^  , . .  are  butputi  of  interest,  a|  influence! 
w'l  and  hence  y'^  and  x^,  and  Eq.  31  can  no  longer  be  ignored.  This  illustrates 
the  reason  why  cohvontionai  transfer  function  blocks  must  be  connected  through 
Isolating  amplifiers  or  their  equivalent.  Isolation  means  that  a^  ■  w|  is  aero 
/-and  that  "Eq. '  r  is  valid. 

If  the  P  matrix  Is  written  for  a  single-input,  single -output  electrleal 
network,  where  It  Is  assumed  that  voltage  E^  Is  the  Input  and  the  current  l2  li 
the  output,  Eq.  29  can  be  reduced  to 

a  »  ll  »  p^jEj  ♦  P22E2  • 

A  comparison  of  these  equations  with,  for  example,  Eq.  276  of  refer¬ 
ence  [8]  shows  that  the  p's  are  (with  the  proper  change  of  signs  and  subscripts) 
identical  to  the  ordinary  "y”  parameters  used  to  describe  a  two-port  communi¬ 
cation  network.  In  fact  the  p's  are  the  short-circuit  driving  point  and  transfer 
admittances  for  the  network.  If  1^  becomes  the  Input  and  E2  the  output,  the 
resulting  p's  will  be  the  "z"  parameters,  or  open  circuit  impedances.  Other 
choices  of  input  and  output  variables  will  result  in  the  "g"  and  "h"  parameters. 
The  transfer  matrix  P  is  essentially  only  a  generalization  of  these  circuit 
parameters  for  multi-port  networks.  (The  "ABC D"  or  general  circuit  param¬ 
eters  which  express  E^  and  in  terms  of  E^  and  while  they  are  a  per¬ 
fectly  accurate  mathematical  way  of  describing  a  system,  cannot  be  considered 
a  special  case  of  the  P  matrix  because  the  corresponding  equations  imply  that  the 
twoconjugate  variables  at  one  terminal  are  bmng  considered  simultaneously 
independent.)  "  : 

This  example  also  serves  to  illustrate  the  fact  that  a  general  n-port  bi¬ 
lateral  network  can  be  represented  by  z"  different  transfer  matrices, 
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depending  upon  the  choice  of  dopendent  and  independent  variable*.  The  P  matrix 
therefore  it  determined  not  only  by  the  physical  ayMtem,  but  alio  by  the  way  In 
^which'it'is'to be'uied. 

0.  A  Matrix  Dofinitibh  for  Output  Interaction 

Since  the  outputs  y^  have  been  considered  dependent  variabies.  we  can¬ 
not  discuss  external  forcing  of  those  variables  without  destroying  the  validity  of 
the  essentially  unilateral  transfer  matrix.  However/  If  y|  is  the  position  of  out¬ 
put  shaft  i,  if  we  apply  enough  torque  to  shaft  l.we  can  produce  whatever  dlaturb- 
ance  is  desired  in  y^.  The  effect  of  a  disturbance  at  output  shaft  1  is  really  the 
effect  of  the  secondary  input  conjugate  to  y|.  If  a  signal  w^  does  not  affect  yj/ 
and  a  signal  Wj  does  not  affect  y^,  then  the  outputs  y^  and  yj  are  non-interact¬ 
ing.  (In  the  position  example,  if  y^  is  obliged  to  take  up  a  position  externally, 
the  torque  necessary  to  accomplish  this  does  not  disturb  yj  and  vice  versa.  ) 

The  P-canonical  structure  of  reference  [2]  is  defined  as  a  system  where 
each  output  depends  upon  all  inputs,  but  is  not  affected  by  disturbance  to  any 
other  output.  This  condition  is  satisfied  If  matrix  B  of  the  preceding  section 
is  diagonal.  The  matrix  A  will  be  identical  to  the  matrix  P  of  reference  [2]. 

In  the  V-canonlcal  structure,  every  output  Interacts  with  every  other 
output.  In  the  matrix  representation,  B  will  have  no  non-zero  terms.  This 
condition  means  that  a  disturbance  applied  to  w^,  the  secondary  input  conjugate 
to  y^,  will  alter  the  values  of  all  other  y's,  as  well  as  that  of  y^. 

The  H-canonical  structure  is  defined  in  reference  [4]  for  a  system 
having  more  outputs  than  inputs.  If  m  -  n  =  i,  it  is  assumed  that  y^  .  .  . 
are  not  affected  by  output  disturbances,  but  that  y^  ^  j  affected  by 

a  disturbance  in  any  output.  This  condition  results  in  a  B  matrix  of  the  fol- 
.'.'lowing'.form:  ■■ 
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B 


P2(n  +  1)  P2(n  +  2)°  *  *  ^  ‘  *  *  ’  "  *  *  *  *  I  *  *  ° 


0  ....  .  .  ..;....  .0 


+  1)  *  "  *  +  0 


P(l  4*  IHn  +  1)  •  •  "  P(l  +  l)(n  +  i  +  I)  fO  +  l)(n  +  m) 


Pm-  • 


+  »«)J 


(33) 


Note  that  while  the  V-canonlcal  structure  Is  limited  to  systems  having  the  same 
number  of  outputs  as  inputs,  there  is  no  reason  to  require  the  presence  of  zeroes 
in  matrix  B,  regardless  of  the  relation  between  m  and  n.  The  P  matrix  always 
relates  n  ■¥  m  dependent  variables  to  n  -f  m  independent  variables,  and  it  is  a 
matter  of  physical  and  mathematical  indifference  whether  n  equals  m.  Com* 
pletely  interrelated  systems  (those  whore  a  disturbance  at  any  output  termin.al 
is  reflected  in  all  outputs)  can  certainly  exist,  even  though  the  number  of  out* 
puts  is  greater  than  the  number  of  inputs. 

In  the  case  of  the  P-canonical  structure,  the  matrix  B  gives  exactly  the 
same  information  as  the  assumption  about  structure,  viz. ,  the  outputs  are 
completely  non-interrelated  in  the  sense  considered  here. 

In  the  V-canonical  structure  the  block  diagram  does  not  give^  quantitative 
information  about  the  effect  at  y^  of  a  disturbance  yj.  If,  however,  in  the  P 
matrix  representation  we  let  all  independent  variables  (x's  and  w's)  be  zero 


except  Wj,  then 


But  also 


i(n  + j)  -  j  • 


(34) 


=  Pj(n  >  j)"'j  ■ 
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Therofore,  tho  effoct  oh  oX  •  diiiturbhnco  at  the  Jth  terminal  produc- 
Ing  . n  output  ^  ^ 

The  fact  that  both  and  yj  are  regardod  &■  the  effecta  of  a  elngle  hauae 
Wj  dob  a  not  alter  the  conclualon  that  an  output  y|  cbrroapohda  to  a  diaturbance 


P  matrix  can  be  completely  determined  externally  lor  an  unknown 
ayatem  only  if  all  independent  variablea  can  be  aucceaaively  mantpulated  and 
all  dependent  variablea  aimultaneoualy  obaerved.  Howeveri  the  degree  to 
which  P  must  be  known  dopenda  upon  the  application;  If  aome  ol  the  aecondary 
Outputa  cannot  be  observed,  the  corrcaponding  rows  may  be  omitted  from  the 
P  matrix  Similarly,  if  aome  of  the  aecondary  inputs  are  always  aero  in  a 
given  application,  the  corresponding  columns  may  be  omitted, 
f.  Bilateral  versus  Unilateral  Matrices 

In  Section  III  c,  it  was  stated  that  for  an  electric  network,  the  P  ma¬ 
trix  can  be  derived  from  Kron's  Z*  matrix,  and  in  reference  [7]  it  is  shown 
that  for  a  network  of  the  typo  considered,  Eq.  2  can  always  be  solved  for  the 
unknown  quantities,  whichever  they  may  be.  This  Is  equivalent  to  saying  that 
whatever  variables  aire  taken  to  be  Inputs  and  outputs  (subject  to  the  constraint 
that  the  current  and  voltage  at  one  node  pair  are  not  both  inputs  or  outputs) 
the  corresponding  P  matrix  can  be  found. 

But  the  P  matrix  represents  n  +  m  equations  in  2(n  +  m)  variables;  if 
the  specification  of  any  n  +  m  uniquely  determines  the  other  n  +  m  (and  the 
above  statement  implies  that  it  does  for  electric  networks  or  their  analogs), 
then  the  equations  represented  by  P  must  hold  regardless  of  which  set  of  yaria 
bles  is  considered  to  be  unknown.  In  other  words,  to  go  from  the  ^  matrix  ) 
representing  one  set  of  Inputs  to  the  P  matrix  of  any  other  set,  it  is  only 
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noceipary  to  solve  the  equations  for  the  new  unknowns.  Thu8,although  Eq.  29 
was  written  under  the  assumptldn  that  the  terms  of  V  were  results  (that  Is, 
determined  by  the  independent  variables  U  and  the  system  equitibne),  it  cort- 
tinuos  to  be  valid  even  if  some  of  the  u*s  become  dependent,  and  some  of  the 
v*s  independent.  For  example,  V  ■  P“^U. 

In  a  sense,  then,  unless  the  system  roprosented  by  P  cannot  be  simu¬ 
lated  by  an  electric  network,  the  P  matrix  description  is  bilateral  —  its  equa¬ 
tions  are  valid  regardless  of  the  choice  of  input  and  output  variables.  The  A 
matrix,  being  only  a  part  of  the  P  matrix,  is  not  bilateral.  Y  *  AX  implies 
that  W  *  0,  while  X  *  A*  Y  implies  that  2  »  0.  Hence  A*  Is  not  generally  equal 
to  A"^.  It  is,  at  least  in  many  cases,  this  reason  which  causes  single-varia¬ 
ble  transfer  functions  to  be  limited  to  signals  passing  in  one  direction, 
g.  State  Vector  Method;  Controllability  and  Observability 

The  P  matrix,  while  providing  information  regarding  output  Interaction, 
suffers  from  the  same  disadvantages  as  the  A  matrix  as  far  as  initial  conditions 
are  concerned.  The  state  vector  method, on  the  other  hand, is  p.irtlcularly 
suited  to  problems  involving  initial  conditions  such  as  those  arising  in  optimal 
control  theory.  For  time-varying  systems  where  Fourier  transform  methods 
fail  and  non-linear  systems  where  the  superposition  principle  is  Invalidated, 
the  designer  has  recourse  only  to  the  state  vectors  or  differential  equation 
method  of  representation. 

By  definition,  when  the  system  is  specified  by  equations  of  the  form 

'  '  ,x  =  Ax  +  Bu  .  '-r:;,', 

y  =  Cs.  +  Du 

the  system  is  assumed  to  be  completely  defined  by  the  "n”  state  variables. 

This  implies  that  the  designer  has  access  to>  all  points  within  the  system  or 
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apriori  knowledgo  of  tho  order  of  the  syatom.  If  every  ateto  of  the  ayatom  ia 
Cdhtrbliablo  then  the  atato  Voctor  can  bo  tranaferred  from  a  given  atato  to  any 
other  doalred  atatb.  it  miiat  bo  remembered  that  any  non-aingular  linear  traha 
formation  Z  ■  PX  yiblda  a  different  but  equivalent  aot  of  atato  variabieg 
For  the  timbdnvariant  situation  auch  a  tranaformatioh  does  not  affect  stability 
or  controllability  of  the  syatom.  While  the  state  vector  rbpreaontation  is  a 
.mathematical  model  of  the  syatom  which  describes  it  completely  It  la  hot  unique 
and  the  state  variables  and  consequently  the  equations  themselves  are  chosen 
to  simplify  manipulations. 

Like  the  A  matrix^  the  state  vector  method  is  uhidlrcctional  in  charac- 

■ter'  '■ 


X{t)  a  ♦ 


t)  B(t)  U{t)  dr. 


X(t)  may  be  determined  from  the  initial  condition  X^  and  tho  input  U(t)  to  the 
system.  Every  state  variable  is  a  dependent  variable  and  may  be  affected 
only  by  changing  the  input  U(t).  Hence, output  interaction  interpreted  as  the 
effect  of  varying  one  output  on  any  other  output  cannot  be  described  using  this 
representation.  On  the  other  hand,  the  homogeneous  solution  f  X^  gives  com¬ 
plete  information  regarding  the  effect  of  an  initial  condition  of  state  variable 
x.^  on  the  response  of  the  system  in  the  absence  of  an  input  U(t). 

A  single-variable  system  which  can  be  represented  by  an  nth-order^^^^^^  ^^^^^^^^^^ 
differential  equation  may  also  be  expressed  in  state  vector  notation  in  the  form 


shown  on  the  following  page,  where  Xj  is  the  output  of  the  system  and  •  •  • 
are  its  (n  -  1)  derivatives.  In  the  multivariable  case,  the  state  vector  X(t) 


consists  of  n  corripbnents  of  which  only  q  are  the  actual  outputs  of  the  system 
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whlle  (n  -  q)  are  secondary  state  variabloa,  some  of  which  are  linear  combina¬ 
tions  of  the  derivatives  of  those  outputs.  The  remaining  ones  Boprosent  unob- 
servaiblo  state  variables  which  do  hot  affect  the  terminal  behavior  of  the  system. 
Under  silro  Initial  cohditlohs,  the  state  vector  repreaohtatlon  may  be  convertad 
into  an  (n  X  r)  matrix  of  transfer  functions  by  the  simple  process  of  taking  trans¬ 
forms.  Of  this,  a  (q  X  r)  submatrix  is  the  same  as  the  A  matrix  of  the  system 
When  the  latter  is  Specified  and  the  system  has  no  unobseryable  state  variables, 
the  remaining  (n  •  q)  rows  of  the  (n  x  r)  matrix  are  determined  by  the  arbitrary 
choice  of  the  secondary  state  variables. 

.  AX-fflU 


When  the  system  has  unobservable  state  variables  the  transfer  matrix 
A  contains  no  Information  about  them.  The  actual  order  of  the  system  Cannot 
be  predicted  using  any  transfer  function  represehtatibn.  In  other  words,  when 
the  transfer  matrix  is  used  it  is  tacitly  assumed  that  unobservable  and  uncon- 
trqjllable  state  variables  are  unimportant.  Probl<*m8  in  which  controllability 
and  observability  do  play  an  important  part  must  be  handled  by  the  state  vector 

■'IV.  Conclusion 

The  aim  of  the  foregoing  discussion  has  been  a  clarification  of  the 
limitations  and  advantages  of  various  methods  of  describing  a  multivariable 
system.  The  results  can  be  summarized  in  the  following  statements. 
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(1)  Tho  convchtio^^^^  matrix, .  or  A  matrix,  describes  th<j  ; 

terminal  behavior  of  a  iihear,  time -Invax'lant  system;  this  dcscrlptlbh  is 
valid  for  a  glv*jn  set  of  input  and  output  varlubles.  The  A  rnatrlx  provides 
tio  Information  about  thb  response  of  the  system  when  ah  externaVsigha!  Is 
applied  at  any  bulput.  It  therefore  does  not  define  output  interaction.  Since 
the  irdnsfer  funciibn?  are  obtained  by  dropping  terms  in  the  transformed 
dlffc rentia)  et|ualioh  coiUalnlng  the  effect  of  initial  conditions,  properties 
like  tortrollablllty  4n«l  obscrvahllUv  must  be  considered  irrelevant  for  sys¬ 
tems  represci  ted  only  by  a  transfer  matrix. 

The  A  matrix  is  thus  best  suited  to  problems  In  which  the  stc.idy 
state  hrhavior  or  transient  response  with  /c ro  initial  conditions  la  of  pri¬ 
mary’  import.ance,  and  where  only’  tcrmln.il  variables  are  of  interest. 

(2)  The  dlvisicf:  of  a  system  represented  by  -i  tr-msfer  matrix  Into 
slrgir-viriable  subsystems  described  by  tr.insfer  functions  dors  rot  provide 
any  additicnal  i:\fwrm-nlon  about  tcrml.n.al  beh.avior.  In  pa rtirular,  output 
intcr.»ctlon  cannot  be  defined  by  mean. h  of  such  subdivision.  Exactly  as  In 
the  c. ISC  of  single -v.irl.ih!c  systems,  all  configur.ation.s  of  subsystems  arc 
eqvuvaltnl  mathcm.it ic :il  mudel.s  of  the  physical  sv.stem;  the  choice  of  a 
parlicu’ar  cor.figurution  tlepencis  o.n  convenience  of  m-inipul.atior..  ^ 

A  similar  method  of  defining  a  mathematical  structure  is  thr  selection 
of  a  particular  set  of  state  variables;  an  infinite  number  of  state  vector  rep¬ 
resentations  are  po.s3ibIe  for  any  linear  time-invariant  system  and  the 
choice  of  any  one  of  these  is  again  a  matter  of  convenience. 

(3)  Prjbierns  involving  signal  now  in  two  directions  at  a  given  ter  ¬ 
minal  can  bo  studied  by  moans  of  the  expanded  transfer  matrix,  or  P  matrix. 
Thi.3  matrix,  while  it  shares  the  inability’  of  the  A  matrix  to  take  into  account 
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non-zoro  initial  conditions,  can  bo  usoful  iii  situations  where  output  interac¬ 
tions  or  non-isolated  intercdnheCtions  of  systems  are  important.  ^ 

(4)  Of  the  methods  considered  in  this  report,  the  state  vector  repre¬ 
sentation  is  the  only  way  of  describing  time-varying  or  noh-llnear  multivaria¬ 
ble  8ysterhs>  since  the  fro^  domain  methods  are  no  longer  valid.  This 
approach  provides  cbmplote  information  regarding  the  response  to  an  input 
Vector  U  in  the  presence  of  inltlarconditions.  It  is  most  suitable  for  problems 
in  which  system  behavior  Is  to  be  controlled  directly  In  the  time  domain,  and 
in  the  area  of  optimal  control.  The  use  of  digital  computers  for  obtaining 
approximate  solutions  makes  this  method  especially  attractive.  However, 
this  method  cannot  be  used  for  systems  with  distributed  parameters  because 
the  number  of  state  variables  becomes  Infinite. 

(5)  The  number  of  Initial  conditions  that  have  to  be  specified  to  Obtain 
a  unique  solution  depends  on  the  order  of  the  system,  and  Is  equal  to  the  num¬ 
ber  of  state  variables.  When  the  system  is  described  only  by  a  transfer  ma¬ 
trix,  however,  the  order  of  the  System  is  not  obvious.  Hence,  in  situations 
where  every  state  of  the  system  must  be  independently  controlled,  only  the 
differential  equation  representation  should  be  used.  It  Is  only  In  Such  situations 
that  controllability  and  observability  are  important. 
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